Model Hamiltonian parameters for half-metallic ferromagnets NiMnSb and Cr02 
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Using the recently developed A''th-order muffin-tin-orbital (NMTO) based downfolding technique 
we revisit the electronic properties of half-metallic ferromagnets, the semi-Heusler NiMnSb and 
rutile Cr02. The NMTO Wannier orbitals for the Mn- ci and CT-t2g manifolds are constructed and 
the mechanism of chemical bonding is discussed. The effective hopping Hamiltonian parameters are 
calculated using a NMTO downfolded basis set. We propose model Hamiltonian parameters with 
possibly minimal basis sets for both half-metallic ferromagnetic alloys. 
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I. INTRODUCTION 

Half-metallic ferromagnets (HMF) such as the semi- 
Heusler NiMnSb or the rutile Cr02 are a subject of 
strongly growing interest mainly due to their poten- 
tial applications in spin-dependent electronics, so-called 
"spintoromcs."i2ii4 Concerning their electronic struc- 
ture, half-metallic materials exhibit a gap in one spin 
channel and a normal metallic behavior for the op- 
posite spin bands. This means that electrons at the 
Fermi level are 100% spin polarized. Much experimen- 
tal effort is devoted to the semi-Heusler NiMnSb com- 
pound, although highly spin-polarized carrier injection 
has not yet been achieved. At low temperatures super- 
conducting point contact measurements revealed a less 
then 50% polarization of valence electrons. ^'^ Similar val- 
ues were obtained by spin resolved photoemission,^ and 
more recently experiments using the synchrotron radi- 
ation showed at room temperature a 40% polarization 
at the Fermi level.* This discrepancy in comparison with 
the predicted half-metallic behavior is attributed to inter- 
face and surface effects. However, by a proper engineer- 
ing, half-metalicity can be restored at surface/interfaces.-S. 
The experimental evidence for half-metallicity is stronger 
in Cr02: Andreev reflectionjiS superconducting tunnel- 
ingfii photoemissionjiSi and point-contact magnetoresis- of two parameters: U„ 



ical point of view. 

On the theoretical side, HMF have been strongly sup- 
ported by first-principles calculations, based on density- 
functional theory.* These calculations offer a proper de- 
scription of the ground state properties, and are usually 
performed for zero temperatures. One way to approach 
the finite temperature behavior is based on modeling the 
many-body interactions in the real material, therefore 
the evaluation of the model Hamiltonian parameters is 
required. These parameters constitute the starting point 
for a finite-temperature full many-body microscopic de- 
scription. The local density approximation (LDA) elec- 
tronic structure allows us to evaluate the effective hop- 
ping parameters, to which a Hubbard type interaction 
is added to construct the starting Hamiltonian. In the 
multiorbital case the Hubbard Hamiltonian is described 
by 
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where a is the spin index, m, to' are local orbitals at site i. 
The on-site Coulomb interactions are expressed in terms 



U, U, 



tanceji^ all give values of polarization in the range 85- 
100 %. Keizer et al. have injected new excitement into 
the field of half metals, by reporting the existence of a 
spin triplet supercurrent through the strong ferromagnet 
CrOs.i* 

Future spin electronic devices based on HMF will prob- 
ably be expected to work around and above room temper- 
ature, so one of the essential requirements is that these 
ferromagnets should have quite high Curie temperatures. 
Both NiMnSb (Tc=730 K) and CrOs {Ta^m K) fulfill 
this requirement. A second essential requirement for us- 
ing half-metallic materials in practical devices could be 
the understanding of finite-temperature behavior of spin 
polarization both from the experimental and the theoret- 
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and Jmm' = This model Hamiltonian is defined on 
a basis set of local orbitals and thus the microscopic in- 
teractions are local, involving a small number of elec- 
trons and a small number of orbitals. A numerically ex- 
act solution of this model is achieved by the quantum 
Monte Carlo solver in the framework of the dynamical 
mean field theory (DMFT)iSiii We used the recently 
developed LDA-f DMFT methoda^^i^*i^^i^°i^^i^^ in order 
to investigate finite-temperature many-body effects for 
the practically important spintronic materials such as the 
semi-Heuslers NiMnSb, FeMnSbj^ and the zinc-blende 
CrAs (Ref. 25) and VAs.^^ In these materials, due to their 
half-metallic ferromagnetic band structure, the incoher- 
ent [nonquasiaprticle (NQP)] states play an important 



2 



role. The NQP states were considered theoretically for 
the first time by Edwards and HertaSi in the framework 
of broad-band Hubbard model for itinerant electron fer- 
romagnets. 

For the realistic electronic structure of NiMnSb NQP 
states are situated just above the Fermi level for the 
minority spin channel, having a considerable spectral 
weig hti^ In FeMnSb the spectral weight of NQP states 
is enhanced in comparison with NiMnSb, producing a 
drastic depolarization at the Fermi leveli^ For CrAs in 
the zinc-blende structure)^ the spectral weight of NQP 
states was studied in connection with the substrate lat- 
tice parameter. For large substrate lattice parameters 
the Fermi level is situated close to the middle of the 
minority spin gap and the NQP states are clearly vis- 
ible. However, for smaller substrate lattice parameters 
the NQP contribution is negligible. VAs in a similar zinc- 
blende structure is predicted by LDA/GGA to be a nar- 
row gap semiconductor J2, In addition to the presence of 
NQP states, many-body interactions determine the clo- 
sure of the semiconducting majority spin gap, leading to 
a half-metallic ferromagnetic ground statef2£ Therefore, 
a correct prediction of new spintronic materials should 
take into account finite-temperature many-body correla- 
tion effects, which might play an essential role in depo- 
larization. 

In order to further investigate the nature of the NQP 
state, a realistic model Hamiltonian is required. In 
this paper we will use the recently developed massive 
downfolding schemeSS*^ in order to produce real space 
Hamiltonian parameters for both NiMnSb and Cr02 half- 
metallic compounds. 

The paper is organized as follows. Section Hll describes 
briefly the computational details of the A^th-order niuflin- 
tin orbital (NMTO) method. The corresponding subsec- 
tions lll Al and lll Bl prcscnt the results of downfolding onto 
the Mn-d, Cr-t2g manifolds, respectively, and gives the 
values of the effective hopping parameters. In the case of 
Cr02 we compare the results of the matrix elements of 
the effective hopping Hamiltonian in two distinct Wan- 
nier orbitals basis sets: the one which describes the full 
Cr-t2g manifold and a second in which dxy and dyz±zx 
orbitals could be described individually. The values of 
the effective Coulomb interaction parameters are given 
in Sec. IIIII In the summary we discuss the construction 
of possible model Hamiltonians with a minimal basis set 
for both half-metallic fcrromagnets. 



II. EFFECTIVE HOPPING PARAMETERS 



TABLE I: The muffin-tin radii for NiMnSb and Cr02. The 
radii of the empty spheres E are also given. The second raw 
coresponds to the LMTO basis sets used in the self-consistent 
calculation of LDA potential. (/) means that the i-partial 
waves were downfolded within the LMTO- ASA -I- cc. 
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The truly minimal set of symmetrically orthonormalized 
NMTOs is a set of Wannicr functions. In the construc- 
tion of the NMTO basis set the active channels are forced 
to be localized onto the eigenchannel Him, therefore the 
NMTO basis set is strongly localized. 

Fourier transformation of the orthonormalized NMTO 
Hamiltonian, H^^^(k), yields on-site energies and hop- 
ping integrals. 
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in a Wannier representation, where the NMTO Wannier 
function |xR.m) orthonormal. 

The matrix element between orbitals m! and m, both 
on site R'=R=0, is i^^, and the hopping integral from 
orbital m! on site R'=0 to orbital m on site R=(a:;, y, z) 

j_xyz 

IS ijv ™ . 




In the present paper we use the NAITO method'^"-'^^ for 
generation of localized Wannier functions. The NMTO 
method can be used to generate truly minimal basis sets 
with massive downfolding technique. Downfolding pro- 
duces minimal bands which follows exactly the bands ob- 
tained with the large basis set. In the case of NiMnSb 
and Cr02, Mn-d and Cr-i2g form the minimal basis set. 



FIG. 1: (Color online) Clt structure with the fee Bravais lat- 
tice (space group _F43m). Mn (green) and Sb (purple) atoms 
are located at (0, 0, 0) and (i, |, 1) forming the rocksalt 
structure arrangement. Ni (orange) atom is located in the 
octahedrally coordinated pocket, at one of the cube center 
positions (|, |, |) leaving the other (i, i, j) empty. This 
creates voids in the structure. 
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FIG. 2: (Color online) The band structure of non-spin- 
polarized NiMnSb calculated with the full basis is given in 
thick (green) line. The thin (red) band have been calculated 
with a Mn-d NMTO basis set. Fermi energy Ep is set to be 
zero. The high-symmetry points are VK(i,l,0), I/(i,i,i), 
r(0,0,0), X(0, 1,0), if(f , |,0), in the W-L-T-X-W-K-T line 
and X(0,0, 1), [/(i,i,l), L(i, |,i), in the X-U-L line. En- 
ergy mesh used for downfolded calculation is given in the right 
of the band structure with a unit of eV. 



The LDA potential is generated with the Stuttgart 
TB-LMTO-ASA code (the LMTO-ASA including the 
combined correction) i^SiS*^ NMTO calculations are per- 
formed using the generated LDA potentials. The radii of 
MT potential spheres and the LMTO bases used in the 
calculation for NiMnSb and Cr02 are listed in Table 
For the detail of the calculation, see the appendixes of 
Rcf. i3i. 



Downfolding onto the Mn- ci manifolds in 
NiMnSb 



The intermetallic compound NiMnSb crystallizes in 
the cubic structure of MgAgAs type (Clfc) with the fee 
Bravais lattice (space group F43m — Tj). The crystal 
structure is shown in Fig. ^ This structure can be de- 
scribed as three interpenetrating fee lattices of Ni, Mn, 
and Sb. The Ni and Sb sublattices are shifted relative 
to the Mn sublattice by a quarter of the [111] diagonal 
in opposite directions. In the present calculation the ex- 
perimental lattice constant of NiMnSb (a=5.927 A) is 
used. 

A detailed description of the band structure of semi- 
Heusler alloys was given using electronic structure calcu- 
lations and tight-binding model analysis and 
we briefly summarize the results. The key points, which 
determine the behavior of electrons near the Fermi level 
for the half-metallic property, are the interplay between 
the crystal structure, valence electron count, covalent 
bonding and large exchange splitting of the Mn-c? elec- 
trons. 




rx u Lw 



FIG. 3: (Color online) The NMTO band structure calculated 
with a Mn-d NMTO basis set. The bands have been decorated 
with (a) t2g and (b) characters. About symmetry points 
and Ef, see Fig.|5| 



1. Chemical bonding and the Mn-d Wannier orbitals 

In the nonmagnetic phase the Clb compounds have 
a d-d gap resulting from covalent hybridization of the 
higher-valent transition metal (Ni) with the lower-valent 
transition metal (Mn). For the minority spin gap open- 
ing, not only the Mn-d-Sb-p interactions, but also Mn-d- 
Ni-d interactions have to be taken into account. More- 
over the loss of inversion symmetry produced by Clfc 
structure (the symmetry lowering from Oh in the L2i 
structure to Td in the Clb structure at Mn site) is an es- 
sential additional ingredient. All the above interactions 
combined with the symmetry lead to a nonzero anti- 
crossing of bands and to the gap opening. The existence 
of sp-valent Sb is crucial to provide stability to this com- 
pound. 

Based on the understanding of the bonding in the 
NiMnSb, we propose a downfolding scheme in which all 
orbitals of all atoms except Mn-d is downfolded. 

In Fig. 121 the non-spin-polarized band structure of 
NiMnSb is calculated with the full basis set (thick green 
line). Mn d, Ni d, and Sb p states are lying between 
—0.5-2.5 eV, —3 — 1 eV, and —6.5 — 3 eV, respectively. 
Sb s states are sitting around —12 — 10 eV, that is not 
shown in Fig. |21 There is an excellent agreement with 
the previous calculations The thin red band has been 
calculated with the downfolded basis set which includes 
only Mn-d orbitals. The energy mesh used in the down- 
folded calculation is given to the right of the band struc- 
ture. The two sets of bands are identical. In Fig. |3| 
the bandstructures of NiMnSb with the orbital character 
projected on to NMTO t2g and eg Wannier orbitals are 
shown. The fatness associated with each band is pro- 
portional to the character of the orbital. The strong 
hybridization between t2g and eg states is clearly seen. 
NMTO Mn-d Wannier orbitals are shown in Fig. |3| The 
triply degenerate manganese t2g orbitals are very compli- 
cated due to the hybridization with Ni-d and Sb-p states. 
The dxy orbital at Mn site is deformed by antibonding 
with the Ni-d state directed tetrahedrally to [111], [111], 
[111], and [ill]. The same Ni-d orbitals couple with Sb- 
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FIG. 4: (Color) NMTO Mn-d Wannier orbitals of NiMnSb. Ni is orange, Mn is green, and Sb is purple. Red (blue) indicates a 
positive (negative) sign. Upper panel: t2g orbitals; dxy (left), dyz (middle), dzx (right). The triply degenerate t2g orbitals can 
be obtained by the permutation of axes. Lower panel: eg orbitals; d^z2_i (left), d^2_y2 (middle). These Eg orbitals are doubly 
degenerated. 

p states. The direct Mn-dxy-Sh-p tt coupling is not seen 

since the distance is d(Mn-Sb):d(Ni-Sb)=l:^. There- 
fore the Ni-d-Sb-p interactions are more favorable. The 
dispersion of the Mn t2g bands is mainly due to hopping 
via the tails of Sb-p and Ni-d orbitals. On the other 
hand, the second nearest neighbor (NN) d-d hopping of 
t2g orbital is small. The Cg orbitals at Mn site are much 
easier to understand: they point towards Sb atoms, and 
a strong pda coupling between Sb-p and Mn-Cg states is 
seen. This induces large second NN d-d hoppings. 



2. Effective hopping matrix elements in the basis set of 
Mn-d Wannier orbitals 

In the many body picture the Mn t2g and Cg constitutes 
the active orbitals which are responsible for the low en- 
ergy physics, having fluctuation in occupation and spins. 
The effective hopping Hamiltonian matrix elements built 
up with these active orbitals are as follows. 

NMTO basis set: 

\X^) ^ {\xy),\yz),\zx),\3z'-l),\x'~y')} . (3) 



The on-site term: 
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where the unit is meV, Ep = 0, and hopping inte 
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to the second NN are shown. The on-site term t[ 
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FIG. 5: Full basis set spin-polarized (ferromagnetic) bands 
for NiMnSb; majority spin (left) and minority spin (right). 
About symmetry points and Ef, see Fig.|5| The similarity of 
the spin-polarized majority with the non-spin-polarized bands 
is evidenced. 



diagonal, t2g and Cg orbitals are triply and doubly degen- 
erated, respectively. The crystal-field splitting between 
t2g and eg orbitals is ^ 74 meV. There are 12 first NN 
and 6 second NN hoppings. Only one hopping integral 
at each NN is shown, but all the hopping integrals can 
be derived from proper unitary transformation due to the 
crystal symmetry. For details, see Ref. 's^. The hopping 
between t2g and Cg orbitals is strongly influenced by the 
presence of tails belonging to the Sh-p and Ni-c? orbitals. 
Due to these tails large values of t2g to eg hoppings are 
evidenced in Eqs. © and ©. We mention that these 
hoppings should be small for the case of "pure" t2g and 
eg orbitals. In addition, due to the presence of Sh-p and 
Ni-d tails t2g and eg orbitals are not divided clearly as 
seen in Fig. |31 Therefore they should be treated equally. 
Hoppings further than third NN are small; for instance, 
the largest values are 66 (third) and 56 meV (fourth), 
a small number of matrix elements of the hoppings are 
less than ~30 meV, and the others are almost zero in the 
third and fourth NN hoppings. Further NN hoppings are 
negligible. 

In this procedure, we obtained the nonzero hopping 
matrix elements between the all Mn-d orbitals, in the 
downfolded representation. The numbers of independent 
parameters of the effective model are equal to the number 
of matrix elements from Eqs. Q)-®. 

We now consider the spin-polarized case. Systems with 
18 valence electrons per unit cell are semiconductors, but 
when they contain more than 18 electrons (22 electrons in 
NiMnSb, that is, Ni Sd^is^, Mn SdMs^ and Sb 5s^5p^), 
antibonding states are populated. Therefore, the non- 
magnetic phase becomes unstable and the magnetic state 
can be stabilized. The large exchange splitting of the Mn 
atom (producing a magnetic moment of around ^3.7 (Ib) 
is crucial to induce a half-metallic property. As we saw 
before, for the nonmagnetic case, the t2g and eg states 
of Ni are situated around 2 eV below Ep. In the spin- 
polarized calculation their position is slightly changed, 
therefore the exchange splitting on Ni is not large. The 
actual magnetic moment calculation gives a value around 
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FIG. 6: (Color online) Cr02 (rutile) structure. Crl (green) 
and Cr2 (orange) are located at (0, 0, 0) and (i,^,^). Cr 
atoms are octahedrally coordinated by oxygen atoms (purple). 
The local coordinate system is used for each Cr atom; xi = 
ib - ^c, yi = -la +^h+ ^c, zi - -I- ^b, 

and X2 = -^a - 5b - ^c, y2 = -^a - 5b + -^c, Z2 = 
— -^a-l- -^b. xi,2 and yi,2 are approximately point to O 
atom, and zi,2 are exactly point to O atom. The local axis is 
transformed into each other by a rotation of 90° around the 
crystal c axis. 

The non-spin-polarized result has a striking resem- 
blance to the majority spin-polarized calculations, pre- 
sented in Fig. [S] Kulatov et al. explained half-metallicity 
of NiMnSb and Cr02 by the extended Stoner factor cal- 
culations in the rigid-band approximation; the minority 
spin-band gap opens due to the exchange splitting, which 
shifts minority bands, so they become empty. According 
to our results, the NMTO antibonding orbitals should 
be a good description for the empty Mn-c? states in the 
minority channel. 

All the above results suggest that a minimal basis set 
which captures the essential physics in NiMnSb can be 
constructed using Mn t2g and eg states. Therefore, the 
non-spin-polarized result with Mn-d orbitals can be a 
good starting point for the many-body calculations. 

B. Downfolding onto the Cr-d manifolds in Cr02 

Chromium dioxide Cr02 has a rutile (tetragonal) 
structure with a=4.421 A, c=2.916 A (c/a=0.65958), and 
internal parameter u=0.30534& The Cr atoms form a 
body-center tetragonal lattice and are surrounded by a 
slightly distorted octahedron of oxygen atoms. The space 
group of this compound is nonsymmorphic (P42 /mnm = 
Z?4^). Cr"*"*" has a close shell Ar core and two additional 
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FIG. 7: Full basis set spin-polarized (ferromagnetic) bands 
for Cr02; majority spin (left) and minority spin (right). Ef 
is set to be zero. The high-symmetry points are r(0,0, 0), 
X(0, i,0), M(i, i,0), Z(0,0, i), i?(0, i, i), ^(i, i, i). 



3d electrons. The Cr ions are in the center of the CrOe 
octahedra so the 3d orbitals are split into a t^g triplet 
and an excited eg doublet. With only two 3d electrons 
the eg states are irrelevant and only the t2g orbitals need 
to be considered. The tetragonal symmetry distorts the 
octahedra, which lifts the degeneracy of the tig orbitals 
into a dxy ground state and dy^+zx and dyz-zx excited 
states,**^ where a local coordinate system is used for 
every octahedron (see Fig. A double exchange mech- 
anism for the two electrons per Cr site was proposed-^ in 
which due to the strong Hund's rule corroborated with 
the distortion of CrOg octahedra leads to the localiza- 
tion of the one electron into the dxy orbital, while the 
electrons in the dyz and d^z are itinerant. 

Measurements of the magnetic susceptibility in the 
paramagnetic phase show a Curie- Weiss-like behavior in- 
dicating the presence of local moments j^"* suggesting a 
mechanism of ferromagnetism beyond the standard band 
or Stoner-like model. 

Several recent experiments including photoemission^ 
soft x-ray absorption, resistivity,^^ and optics^ sug- 
gest that electron correlations are essential to the un- 
derlying physical picture in Cr02. Schwara^ first pre- 
dicted the half-metallic band structure with a spin mo- 
ment of 2fjiB per formula unit for Cr02. Later on Lewia^i 
used the plane-wave potential method and investigated 
the energy bands and the transport properties, charac- 
terizing Cr02 as a "bad metal," terminology applied to 
high temperature superconductors, or to the other tran- 
sition metal oxides, even ferromagnets such as SrRuOa. 
A decade later the LSDA-|-f/ calculation^ explained the 
conductivity in the presence of large on-site Coulomb in- 
teractions, and concluded that Cr02 is a negative charge 
transfer gap material which leads to self-doping. Con- 
trary to the on-site strong correlation description, trans- 
port and optical properties obtained within the FLAPW 
method (LSDA and GGA)j5S suggest that the electron- 
magnon scattering is responsible for the renormalization 
of the one-electron bands. In addition, more recent model 
calculations proposed even orbital correlations!^ 



1. Chemical bonding and the Cr-t2g Wannier orbitals 

Chemical bonding in rutile-type compounds including 
Cr02 was analyzed by Sorantin and Schwarz.-- Let us 
summarize their results: one can see that around the 
Fermi level, the bands are primarily chromium 3d states 
of tig manifold, with Cg bands situated higher in energy 
due to the crystal-field splitting. In the spin polarized 
case, the exchange splitting shifts the minority spin-d 
bands above the Fermi level, as seen in Fig. For the 
majority tig bands the Fermi level is lying in a pseudo- 
gap. Oxygen p-chromium d, hybridization creates both 
bonding and antibonding hybrid orbitals, with the bond- 
ing orbital appearing in the occupied part and the anti- 
bonding hybrid orbital remaining in the Cr tig manifold. 
Half of the dyz and dzx components of tig are pushed up- 
ward, by antibonding, that explains the dominance of dxy 
character in the spin density. The non-magnetic density 
of states shows a sharp peak at the Fermi level, which 
signals the magnetic instability according to the usual 
Stoner argument. 

By the procedure of NMTOs we can easily gain a more 
complete picture about the chemical bonding. In Fig. |H| 
the non-spin-polarized band structure of Cr02 is calcu- 
lated with the full basis set (thick green line). Cr Cg, Cr 
tig, and O p states are lying between 1.5-5 eV, —1-1.5 
eV, and —7.5 — 1.5 eV, respectively. There is an excel- 
lent agreement with the previous calculationsi^2ii2iSiS 
The thin red bands have been calculated with the down- 
folded basis set which includes only Cr-d orbitals. The 
energy mesh used in the downfolded calculation is given 
to the right of the band structure. NMTO Cr-tig Wan- 
nier orbitals at Crl (0, 0, 0) and Cr2 (i, i, i) are shown 
in Fig. 13 All tig orbitals form antibonding pdi: cou- 
pling with 0-2p states. These antibonding NMTO Wan- 
nier orbitals capture the essentials to describe the half- 
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FIG. 8: (Color online) The band structure of non-spin- 
polarized Cr02 calculated with the full basis is given in thick 
(green) line. The thin (red) band have been calculated with 
a CT-t2g NMTO basis set. Energy mesh used for downfolded 
calculation is given in the right of the band structure with a 
unit of eV. About symmetry points and Ep, see Fig.Q 
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FIG. 9: (Color) NMTO Cr-t2g Wannier orbitals of Cr02. Crl is orange, Cr2 is green, and O is purple. The local coordinate 
systems are used for each atom. Red (blue) indicates a positive (negative) sign. Upper panel: Crl t2g orbitals. Lower panel: 
Cr2 t2g orbitals. d^cy (left). This orbital is singly degenerated, dy^-zx (middle) and dy^+zx (right). These two orbitals are 
nearly degenerated. 



metallicity of Cr02. 

Note that the local coordinate system is usediS For in- 
stance, the z axis at Crl site points along the [110], while 
the z axis at Cr2 sites points towards the [110] direction. 
In the following we describe the essential features of the 
Cr-c? and 0-2p orbital couplings in the t2g orbitals on Crl 
site: 

(i) CrsO cluster. The oxygen atoms yield a CraO clus- 
ter with three surrounding Cr atoms. The Crl atom has 
six CraO clusters with two different types; for instance, 
one CraO cluster is formed by Crl at (0,0,0) and (0,0,1), 
Cr2 at (-i,i,i), and O at (u - \,\ - u,\) in the (110) 
plane. Another type of CrsO cluster is formed by Crl at 
(0,0,0), Cr2 at (i,i,±i), and O at (m,u,0) in the (IlO) 
plane. The O is sitting in the center of a triangle formed 
by the three coplanar Cr nearest atoms. The CraO clus- 
ter is indicated by thick cyan lines in Fig.^ In the CraO 
unit, one can see one antibonding pdi: coupling between 
Cv-t2g and 0-p states and two bonding pda couplings be- 
tween Cr-Cg and 0-p states. On one hand, the Cr-e^ and 
Cr-t2g orbitals lying within the CraO plane couple to the 
in-plane 0-p orbital. On the other hand, the out-of-plane 
O p orbital component, perpendicular to the CraO plane, 
cannot couple to the Cr-e^ state due to orthogonality. 



(ii) dxy orbitals. The Crl atom experiences a bond- 
ing coupling between its dxy orbital and the Cg {d^z'^-i) 
orbital located on the other six nearest Cr atoms on the 
(110) plane via CraO cluster. The bonding is realized via 
the tails of 0-2p {p^- and py-like) orbitals. eg tails be- 
longing to the first NN Crl (0,0, ±1) and second NN Cr2 
(— ^,^,±i), (i,— i,±i) atoms are visible in the (110) 
plane, where the dxy orbital is situated. Further, the 
dsz^-i state belonging to the Crl (0,0,1) atom forms 
a pdcr-type bonding with both O Px orbital located at 
(m — — u,^) and with O Py situated at {^ — u,u—^,^). 
The dxy orbital, on the contrary, forms an antibonding 
coupling with these O px and O py states. Cr2 atom sit- 
uated at (—5,5,5) involves its d^.z'^-i orbital into a pda 
bonding with O px orbital situated at (u — — ii,i). 

(iii) dyz-zx orbitals. The O atoms situated at (u,m,0) 
and (— u,— u,0) intermediate a pda bonding, via their 
Px-y orbital, between the dyz-zx orbital at Crl (0,0,0) 
and the eg orbitals at Cr2 (5, 5, ±5) and (—5,— 5, ±5) 
atoms. This bonding orbital is situated in the (110) 
plane similarly to the dyz-zx orbital. However, its bond- 
ing strength seems to be weaker than the ones of the 
dxy orbitals, due to a larger distance between O and Cr2 
atoms. The O Pz orbitals belonging to the atoms situated 
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at 



i,±i) and {u 



1 1 

2 ' 2 



M,±i), form a pdn 



coupling with dy^-zx orbital. Therefore the O Pz orbital 
perpendicular to the plane of the CraO cluster does not 
overlap with the Cg tails of Crl (0,0, ±1) atoms sitting 
along the [001] direction. 



(iv) d. 



yz+zx 



orbitals. The Cg tails can not contribute 
at all to the dyz+zx orbital. Because O pz orbitals at 
(i — u,u — ^,±5), {u — i,i — u,±i) and O Px+y orbitals 
at (u,M,0), {—u,—ufi) are situated perpendicular to CraO 
clusters, the coupling between O p (jpz and Px+y) states 
and Cg states at surrounding Cr atoms is not allowed due 
to the orthogonality. 

We discussed in the above points (i)-(iv), the direct or 
mediated interactions between Cr-d and 0-p or between 
the Cr t2g and Cg states. As it is already known ^^i^^i^SiS 
there is a significant difference between the t2g and Cg 
orbitals, however the above analysis in the framework 
of NMTO technique shows that their interplay consti- 
tutes an important ingredient not only for the crystal- 
field splitting of t2g states, but also for the general bond- 
ing in the rutile structure. The t2g orbitals form the basis 
set in which the effective hopping Hamiltonian matrix el- 
ements are evaluated. These results are presented below. 



Effective hopping matrix elements in the Cx-t2g Wannier 
orbitals 



The hopping integrals t with t2g Wannier representa- 
tion up to the second NN are as follows. 
NMTO basis set: 



Ix"^) = {\xy), \yz-zx), \yz^'zx)} . 
The on-site term: 



000 



104 
323 
355 



(7) 



(8) 



The first nearest neighbor: 



OOl(Crl-^Crl) 



The second nearest neighbor: 




(9) 



1 i i 
^222 



(Crl^Cr2) 



-4 
32 
-204 142 



where the unit is meV and a local coordinate system for 
each atom Crl and Cr2 in the unit cell is used as seen 
in Fig. El The on-site term and the first NN hop- 

pings are diagonal in a {dxy, dyz-zx, dyz+zx} representa- 
tion. The splitting between these orbitals is due partly 
to the orthorhombic distortion of CrOg octahedra and 
partly to the bonding with the Cg states of nearest Cr 
atoms in rutile structure. As a consequence mainly two 




FIG. 10: (Color online) NMTO downfolded d^y (left) and 
dyz±zx (right) bands for Cr02. Energy mesh used for each 
calculation is given in the right of the band structure with a 
unit of eV. For symmetry points and Ef, see Fig. |7| 



crystal-field levels are formed, a single dxy and two nearly 
degenerate dyz±zx orbitals situated with 235 meV at 
higher energies, that can be derived from the difference 
between on-site levels of dxy and dyz±zx orbitals in Eq. 
(jHl). A possible ingredient in determining the position 
of the nearly degenerate dyz±zx orbitals with respect to 
the single d^y orbital is the bond length between Cr and 
O atoms. In the calculations presented in this paper a 
longer bond length along z direction (c?_l — 1.91 A) is 
used in comparison with the bond formed in the xy plane 
having a value of (dy = 1.89 A). In another experiment, 
however, the opposite situation is reported, d±_ = 1.89 
A is smaller than d\\ = 1.91 Ai^ In both cases, how- 
ever, the band structure and density of states are almost 
identical. Moreover other oxides (Ti02, VO2 etc.) with 
rutile structure have a similar band structure and density 
of states for t2g states. This fact suggests that the bond- 
ing with Cg state is essential for the crystal-field split- 
ting of t2g states. Similar mechanism for the crystal-field 
splitting of t2g state is realized in V2O3 (Ref. 55) and 
NaCoOa-^ 

There are two first NN between Crl (Cr2) and eight 
second NN hoppings between Crl and Cr2. Only one 
hopping integral at each NN is shown here. The hopping 
integrals are strongly affected by Wannier function's tails 
which are mainly due to downfolded O 2p and Cr Cg or- 
bitals. C?r 



'■xy 
111 

integrals t^ ^ ^ 



xy^xy 

This means that d. 



orbital has a small second nearest hopping 

4 meV although O p tails are present, 
orbital does not have the second 



xy 



(^20) NN hopping path, indicating a very narrow d^y band. 



On contrary, large hopping integrals of the second NN are 
shown in dyz±zx block, dyz+zx orbitals can hop to dyz+zx 
of all direction (|iy"+i^'^y2+23;| = 142 meV), but dyz~zx can 



yz-zx {t1^'^f^y^_^^^{) meV exactly, as a conse- 



not to d. 

quence of orthogonality.) The hopping between dyz-zx 
and dyz+zx (|i|=204 meV) is allowed for a particular di- 
rection, that is, the second NN hopping from Crl dyz+zx 
to Cr2 dyz-zx sitting in (110) plane or from Crl dyz-zx 
to Cr2 dyz+zx sitting in (110) plane. The large values of 
second NN hopping in dyz±zx block are caused by O p 
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tails which are not orthogonal to each other, and induce 
both the itinerant property and the wide band of dyz±zx 
states. dyz±zx orbitals have the second NN hopping path 
via O p tails. Moreover the large ofT-diagonal hopping 
integral of second NN and the small crystal-field split- 
ting between dyz±zx indicate that dyz±zx states should 
be treated together. The Crl dyz-zx orbital is weakly 
coupled with Cr2 dxy {\t\=32 meV) along particular di- 
rections, within the (110) plane. Hoppings towards the 
third NN are small; for instance, the largest values are 
45 meV (third) and 30 meV (fourth and fifth), and the 
others are almost zero. Further away NN hoppings are 
much smaller. 

Exceptions are, however, the eg tails of Wannicr func- 
tion which gives a few large hopping integrals, for in- 

, ,002(Cil^Crl) 

stance, ixy,xy 



-82 meV, t 



meV, and t 



lli(Crl-+Crl) 



llO(Crl^Crl) 
xy,xy 



matrix elements in LDA Hamiltonian between dxy and 
dyz±zx block, that is there is no LDA interaction between 
them, by virtue of our construction of NMTO Wannier 
functions. 

The possibly minimal set of Wannicr function may pro- 
vide new insights. At the first glance one can separate 
the narrow dxy orbitals from the extended dyzztzx 

states. 

The former would be treated as corelike (classical) spin 
S — 1/2, meanwhile in the later dispersive dyz±zx bands 
the Coulomb repulsion would be treated in a usual quan- 
tum many-body way. Such a Kondo-lattice type model, 
for Cr02 would be described by the following Hamilto- 



= -95 



-yz-zx,yz-zx = "50 mcV. Thcsc large hop- 
pings can be seen in Fig. |3 



3. Effective hopping matrix elements in the independent 



1 

2 



(15) 



d^y and d, 



yz^zx 



Wannier orbitals 



As seen in Eqs. (|Sll- HlU|l . the coupling between dxy 
and dyz±zx seems to be weak. The NMTO method may 
be able to pick up such bands independently. The dxy 
and dyz±zx bands are shown in Fig. llUI The Wannier or- 
bitals are more extended due to additional downfolding. 
This affects the hopping integrals as well. The hopping 
integrals t with dxy and dyz±zx Wannier representation 
up to the second NN are as follows. 

NMTO basis set: 

= {{\xy)}, {\yz-zx), \yz + zx)}}. (11) 
The on-site term: 



-24 



poo 



312 
373 



(12) 



The first nearest neighbor: 



rOOl(Crl^Crl)^ 
va'ra 



The second nearest neighbor: 



^i|i(Crl->Cr2)^ 



,002(Crl^Crl 
f'xy^xy 



and t. 



lli(Crl-+Crl) 
yz~zx^yz — zx 



The first term denotes the hopping of the conduction 
electrons between the NN sites i,j with a hoping matrix 
element t™-™ between the dyz±zx orbitals, described by 
the m,m' indices. The dyz±zx electrons and the local- 
ized dxy spins interaction is given by the second term as 
an exchange coupling (J < for antiferomagnetic and 
J > for ferromagnetic coupling). The ferromagnetic 
coupling between the core spins and conduction elec- 
trons, favors ferromagnetic ordering, because the hop- 
ping amplitudes of the conduction electrons reach the 
maximum possible values if the core spins are aligned. 
Many of the Mn-type collosal magnetoresistance materi- 
als (CMR) can be described by such a model. Recently 
a realistic LDA-I-DMFT calculation for such a model is 
carried outiS The d shell of the Mn'^+ in the undoped an- 
tiferromagnetic insulator contains three electrons in the 
t2g orbitals forming a core spin of magnitude S = 3/2, 
which due to strong Hund's rule coupling couples ferro- 
magnetically to one additional electron in one of the eg 
orbitals. For a certain doping CMR materials are fer- 
romagnetic metals because of additional holes in the eg 
conduction bands. 

In contrast to the manganites, in Cr02, the dxy spin 
has a smaller in value S — 1/2, and the splitting of the lo- 
calized dxy orbital with repespect to the itinerant dyz±zx 
one is smaller [« 370 meV, see. Eq. (|12f) ]. 

It was shown recently^^ that for a proper definition 
of the many-body green function in Cr02, the complete 
t2g manifold is required. Moreover a better description of 
polarization can be obtained considering the competition 
between quasiparticle description around the Fermi level, 
and a local moment behavior at higher energies, above 
the Fermi level 

According to our results, NMTO can give a minimal 
model with a possibly minimal set of distinct Wannier 
-39 meV. In addition, there is no dxy and dyz±zx functions. This result does not exclude 



-99 



-26 



-187 
197 




-202 135 



(13) 



(14) 



where the same definitions of Eqs. l(7|l- l(T(Hl are used. 
Since the Wannier orbitals are more extended the hop- 
pings are modified slightly, but not changed signifi- 
cantly. A few large hoppings of far NN are also reduced; 

21 meV, tiy^iy'^^^'^'^^' = -16 meV, 
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the possibility of integrating out the degrees of freedom 
connected to the narrow/extended d^y/ dyz±zx orbitals. 
Physically this would correspond to the dualistic charac- 
ter of the electronic strcutre of Cr02 around the Fermi 
leveliS 



III. EFFECTIVE COULOMB REPULSION 

The other essential component for a model Hamilto- 
nian describing correlations, is the average Coulomb in- 
teraction parameter U . This term act on the diagonal 
part of the effective Hamiltonian and corresponds to the 
screened electron-electron repulsion. Aryasetiawan et al. 
pointed out recently^ that a rigorous way to define this 
quantity can be formulated in terms of path integrals by 
performing a partial trace, over the degrees of freedom, 
that one wants to eliminate. However, in practice for re- 
alistic materials, the elimination of degrees of freedom is 
a very difhcult procedure. 

In order to evaluate the average Coulomb inter- 
action on the d atoms and the corresponding ex- 
change interactions we start with the constrained LDA 
j^gl-jjojj 61,62,63,64,65^ j-^ |.jjjg approach the Hubbard U 

is calculated from the total energy variation with respect 
to the occupation number of the localized orbitals. In 
such a scheme the metallic screening is rather inefficient 
for 3d transition metalSfS^ The perfect metallic screening 
will lead to a smaller value of U. Unfortunately, there are 
no reliable schemes to calculate U in metalSfS^^ therefore 
in our previous works2i2iiS42& we choose some interme- 
diate values of U from 2 to 4.8 eV and J = 0.9 eV. 

In the case of NiMnSb the constrained LDA calculation 
indicates that the average Coulomb interaction between 
the Mn 3c? electrons is about U = 4.8 eV with an ex- 
change interaction energy about J = 0.9 eV as seen in 
Tableim A reduced value of U can be obtained if one con- 
sider that the Mn eg orbitals participate in the screening 
of Mn t2gr - reducing the average Coulomb interaction 
on the Mn atoms to a value of [/ = 4.2 eV. However, our 
DMFT resultSiSi showed that the many-body effects are 
equally essential for t2g and Cg orbitals, therefore a model 
with only Mn t2g, even it may capture the main physical 
results, would suffer of incompleteness. Note that physi- 
cal results for NiMnSb are not very sensitive to the value 
of U, as it was demonstratedi^^ 



TABLE II: The constrained LDA values of the average 
Coulomb and exchange interactions. The second raw core- 
sponds to the results when the Cg orbitals screen the t2g ones. 
This type of screening would be more appropriate for the 
Cr02 case. 





NiMnSb 


CrOa 




U (eV) J (eV) 


U (eV) J (eV) 


t2g and eg 


4.80 0.93 


3.50 0.90 


t'Zg 


4.25 0.93 


3.00 0.87 



Concerning the local Coulomb interaction U , in Cr02, 
we saw from the analysis of the Cr-3ci manifold that the 
higher energy Cg bands makes no noticeable contribution 
to the Fermi level, however they could participate in the 
screening of the t2g orbitalSjSiSS giving the following val- 
ues for [/=3 eV and J=0.87 eV. 



IV. SUMMARY 

Spintronics requires the search for new materials 
such as half-metallic ferromagnets, whose properties are 
commanded by their electronic structure and electron- 
electron correlations. To understand the collective effects 
a clear picture of the interplay between the microscopic 
interactions is necessary. Model Hamiltonians allows a 
controlled reduction of the ah initio information into a 
few dominant material specific hopping and interaction 
parameters. In this paper we have applied the downfold- 
ing procedure within the iVth order mufBn-tin orbital 
(NMTO) method to obtain model Hamiltonian parame- 
ters for two ferromagnetic half metals NiMnSb and Cr02. 

For NiMnSb, the present Wannier orbital results con- 
firm the previous conclusions of Ref. ^ based on the 
gamma point analysis. However based on the whole Bril- 
louin zone analysis we found that d-d hybridization be- 
tween the transition metal atoms is essential for the gap 
formation as the p-d hybridization previously discussedi^ 
Due to the significant Mn exchange splitting, a different 
hybridization takes place in the spin-up and spin-down 
channels, which determines the appearance of the gap 
only for spin-down electrons meanwhile for spin-up states 
a metallic character is evidenced. 

A downfolding calculations which would include the 
Sb-p and Ni-d orbitals, would increase considerably the 
computational demand. However, the physical picture 
would not be very much changed, because an antibond- 
ing Mn d Wannier orbital including Sb p and Ni d char- 
acter as a tail should play an important role for half- 
metallic nature around Fermi level. Therefore, the effec- 
tive Hamiltonian described by Eqs. Q-®, represents a 
good starting point to investigate many-body effects. 

Accordingly many-body correlation effects would be 
of primary importance for the Mn orbitals the Ni one 
could be neglected. In our previous many-body re- 
sults^^i^ we made use of the above conclusions, so cor- 
relation effects were considered only for Mn-d orbitals. 
The LDA-I-DMFT results showed the existence of non- 
quasiparticle states (NQP) stateSLSS in the minority spin 
channel. These states are associated with finite temper- 
ature spin fluctuation processes which play an impor- 
tant role in depolarization. In the majority spin chan- 
nel, FeMnSb (Ref. ^ has a larger DOS at the Fermi 
level in comparison with NiMnSbi^ This contributes to 
a stronger finite temperature depolarization in FeMnSb, 
which increase with the on-site Coulomb interaction. 

For Cr02 the bands around the Fermi level are primar- 
ily chromium id states of t2g manifold, with Cg bands 
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situated higher in energy by the crystal-field splitting. 
The orbitals are further splitted into a single d^y and 
nearly degenerate dyz±zx bands due to the orthorhombi- 
cal distortion of CrOe octahedra. We discussed the direct 
and the mediated Cr-d, 0-p, or Cr t2g-eg interactions. 
Despite the differences between Cr t2g-eg orbitals their 
interaction play an important role not only in character- 
izing the crystal-field splitting, but also in the general 
picture of bonding in the rutile structure. 

The basis set of the Wannier orbitals for Cr02 is 
formed in t2g manifold. In the framework of the NMTO 
technique two distinct downfolding procedures are possi- 
ble. 

In the first approach a full t2g construction is possible, 
i.e., a single d^y and a nearly degenerate dyz±zx orbitals 
represents the basis set. We believe that such a basis set 
can be used in connection with the multiorbital Hubbard 
Hamiltonian described in Eq. (^). 

A second approach is made possible due to a weak d^y 
and dyz±zx coupling. In this case the lowest dxy and the 
higher dyz±zx bands can be derived separately as shown 
in Fig. ^1 For such a situation, downfolded NMTO 



Hamiltonian with a set of those Wannier orbitals {dxy} 
and {dyz+zx, dyz-zx}, is block diagonal, that is, there is 
no hopping between the lowest dxy and the higher dyz±zx 
orbitals. This fact gives a strong impression for a Kondo- 
lattice type model Hamiltonian described by Eq. H15() . In 
a future work we report on the qualitative/quantitative 
differences of such models applied to Cr02. 

Similarly to our previous DMFT results on semi- 
Heuslers22iSi and zinc-blende structures further 
work will include the many body effects in Cr02 and 
the analysis of nature of the NQP states in these classes 
of half-metallic ferromagnets. 
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